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ABSTRACT 

We present Hubble Space Telescope (HST) WFC3 imaging and grism spectroscopy observations of 
the Herschel- selected gravitationally-lensed starburst galaxy HATLASJ1429-0028. The lensing system 
consists of an edge-on foreground disk galaxy at z = 0.218 with a nearly complete Einstein ring of 
the infrared luminous galaxy at z = 1.027. The WFC3 spectroscopy with G102 and G141 grisms, 
covering the wavelength range of 0.8 to 1.7 /im, resulted in detections of Ha+[NII], H/?, [SII], and 
[OIII] for the background galaxy from which we measure line fluxes and ratios. The Balmer line ratio 
Ha/H/3 of 7.5 =b 4.4, when corrected for [Nil], results in an extinction for the starburst galaxy of 
E(B — V) = 0.8 zb 0.5. The Ha based star-formation rate, when corrected for extinction, is 60 zb 50 
M 0 yr _1 , lower than the instantaneous star-formation rate of 390 zb 90 M 0 yr _1 from the total IR 
luminosity. We also compare the nebular line ratios of HATLASJ1429-0028 with other star-forming 
and sub-mm bright galaxies. The nebular line ratios are consistent with an intrinsic ultra-luminous 
infrared galaxy with no evidence for excitation by an active galactic nucleus (AGN). We estimate the 
metallicity, 12 + log(0/H ), of HATLASJ1429-0028 to be 8.49 zb 0.16. Such a low value is below the 
average relations for stellar mass vs. metallicity of galaxies at z ^ 1 for a galaxy with stellar mass of 
^ 2 x 10 11 M 0 . The combination of high stellar mass, lack of AGN indicators, low metallicity, and 
the high star-formation rate of HATLASJ1429-0028 suggest that this galaxy is currently undergoing 
a rapid formation. 

Subject headings: cosmology: observations — submillimeter: galaxies — infrared: galaxies — galaxies: 
evolution 


1. INTRODUCTION 

Dusty star-bursting galaxies, especially those that are 
identified at far-IR/sub-mm wavelengths, have infrared 
luminosities Lir ^ 10 12 — 10 13 L 0 , implying star- 
formation r ates (SFRs) in excess of 200 M 0 yr -1 (see 
review by (jCasev. Narayanan, fc Cooravl 120141 )). As a 
primary contributor to the cosmic far-IR background, a 
significant fraction of cosmic star formation and metal 
production could have occurred in these star-bursting 
galaxies. Due to deep and wide surveys with the Her¬ 
schel Space Observatory (Pi lbratt et al.l 120 10). we now 
have large samples of dusty, star-burst galaxies at z > 1. 
Despite large number statistics our knowledge on the 
physical processes within such galaxies is still limited. 

Traditional studies at optical and IR wavelengths in¬ 
volving nebular lines to probe the interstellar medium 
(ISM) of these dusty starbursts are challenging due to 
high dust extinction. One way to overcome this limita¬ 
tion is to make use of the flux magnification provided 
by gravitational lensing. Sub-mm surveys provide an 


efficient way to select lensed high-redshift galaxies due 
to the negative K-correction of the thermal dust spec¬ 
tral energy distribution (SED) an d the steep faint-end 
slope of the sub-mm source counts (Blain 1996). T he two 
large area sur veys, Herschel -ATL AS (jEales et al.ll20loh 
and HerMES ([Oliver et al.l l2Q12h . have r esulted in suf¬ 
ficiently large samples of lensed galaxies dNegrello et al.l 
120101 : 1 Wardlow et al.l 120131 : iBussmann et alJ l2Q13h from 
which we can find interesting targets for detailed follow¬ 
up observations. 

Here we present results on the rest-frame optical spec¬ 
troscopy of a lensed starburst galaxy to study its nebular 
line emission and line ratios. We make use of the Hub¬ 
ble Space Telescope Wide Field Camera 3 (HST/WFC3) 
grisms for these observations. To detect both Ha and 
H/? over the wavelength covered by WFC3 grisms we 
require lensed galaxies to be at z < 1.6. One feasi¬ 
ble target for WFC3 grism observations from currently 
known l ensed Herschel sou rces is HATLASJ142935.3- 
002836 (IMessias et al.l FUTl (H1429-0028; G15v2.19 in 
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ICalanog et al.l ([2014l) h The lensed galaxy was detected 
in th e GAMA-15 field of Herschel -ATLAS (jEales et all 
2010 ) with Si60/xm=l-l A 0.1 Jy. The lensing models of 
the system using KeckII/NIRC2 laser guide star adaptive 
optics image a nd high-resolution AL MA imaging data 
are pr esented in ICalanog et al.l (j2014f ) and IMessias et al.l 
( 12011 . The system includes a foreground edge-on disk 
galaxy (z = 0.218) with a near complete Einste in ring 
(Fig. 1 ). The lens model in IMessias et al.l ( 2014 ) shows 
that H1429-0028 is comprised of two components with 
a mass ratio of (1 : 2 . 8 t}' 5 ). The t wo components have 
been used to suggest H1429-0028 may be undergoing a 
galaxy merger, but the two components of H1429-0028 
are found to lie on top of each other. This also leaves the 
possibility that the compact bright component of H1429- 
0028 is a starbursting clump or a region within a galaxy. 
The full extent of the galaxy is traced by the extended 
component that is gravitationally lensed to an Einstein 
ring (Fig. 1). The background source is at z = 1.027 
with a total magnifi cation factor o f ti = 7.9 db 0.8 at in¬ 
frared wavelengths (jMessias et al.l 120141 ) . The K-band 
(AB) magnitude of H1429-0028 is 18.2 (jCalanog et al.l 
120141 ) and is at the level that allows grism observations 
with just one or two HST orbits. 

Here we report HST/WFC3 grism spectroscopic ob¬ 
servations of H1429-0028 making use of G 102 and G141 
grism filters, covering the wavelength range of 0.8 to 1.7 
/im. At the redshift of H1429-0028 these observations 
probe the rest-frame wavelength range of 0.4-0 .8 fi m al¬ 
lowing us to measure several photoionization emission 
lines. We use these measurements to explore the proper¬ 
ties of this system in terms of several emission line diag¬ 
nostic diagrams. We also establish the gas-phase metal- 
licity in a star-forming galaxy. The Paper is organized as 
following: In Section 2 we describe the observations and 
our data reduction procedure. In Section 3 we present 
our results related to emission lines and emission line flux 
ratios and discuss them in the context of existing stud¬ 
ies in the literature. We conclude with a summary in 
Section 4. When calculating luminosities we make use of 
the standard flat-ACDM cosmological model with Ho = 
70 km s -1 Mpc -1 and Da=0.73. 

2. OBSERVATIONS 

HST/WFC3 observations of H1429-0028 were com¬ 
pleted with two orbits under GO program 13399 in Cycle 
21 (PLCooray). We obtained a total of five exposures, 
including two direct images and three grism observa¬ 
tions in two filters. The two direct images made use 
of WFC3/F160W and F105W filters for a total of 250 
and 350 seconds, respectively. We obtained G 102 and 
G141 grism observations over 1800 and 2900 seconds, re¬ 
spectively. The G141 grism covers 1.0 to 1.8 /im, while 
G 102 grism covers 0.7 to 1.2 /im. At z = 1.027 these 
observations then allow important emission line studies 
of H1429-0028 involving Ha at 1.33 /im, H/? at 0.985 /im, 
[OIII] at 1.015 & 1.005 /im, and [SII] at 1.364 & 1.361 
/im. Due to the low spectral resolution of order 80 A the 
data do not resolve the [SII] doublet or Ha from [Nil]. 

We made use of the calibrated HST imaging and grism 
data from the CALWF3 reduction pipeline, as provided 
by the Space Telescope Science Institute. The spectra for 
individual objects in t he image were extrac ted with the 
aXe software package (|Kiimmel et al.l 2009j. Briefly, we 



Fig. 1.— Top Left: The three color image of the grav¬ 
itationally lensed system HATLASJ1429-0028 (also G15v2.19 
in ICalanog et all (l20Ti b using WFC3/F105W (blue), F160W 
(green), and KeckII/NIRC2-LGS K s (red) imaging data. Top 
Right: The foregr ound lensing galaxy was modeled using GALFIT 
(IPeng et al.H2002l) and then subtracted from the K s -band imaging 
data. We l abel the bright knots following the scheme that was 
presented in IMessias et all ( 20141 "). Bottom Left: Source plane re¬ 
construction showing the two components of HATLASJ1429-0028 
and the caustic curves. Bottom Right: Image plane reconstruction 
showing the lens model and the critical curve. Note that the bright 
features A, B, C and D are from the bright compact source near 
the inner cusp caustic while the diffuse ring is due to the extended 
source. 
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Fig. 2.— Left: The direct image in each of the WFC3 imaging 
filters oriented along the dispersion direction of the grism. Right: 
The two-dimensional grism images of H1429-0028. The top panel 
shows G141 and F160W images while the bottom panel shows G102 
and F105W images. The upper trace contains the signal from knots 
A+B while the lower trace contains the signal from knot C as well 
as from the foreground lens. 

created an object catalog making use of the broad-band 
F160 W and F1Q5W images w ith the S Extractor pack¬ 
age (jBertin fe Arnoutsl Il996 ). A sky background sub¬ 
traction was performed on the direct as well as the grism 
images. The core aXe marks spectral regions for each 
object in the SExtractor catalog, estimates contamina¬ 
tion from nearby sources, and flat fields each of those re¬ 
gions or beams. A two-dimensional stamp of each grism 
beam is generated and then combined together with mul¬ 
tiple observations of the same stamp to create a final two- 
dimensional image for scientific analysis. The data prod¬ 
ucts include the two-dimensional combined grism stamp 
for each object as well as flux-calibrated one-dimensional 
spectra, contamination estimates, and error estimates. 

We identified emission lines in the one-dimensional 
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Fig. 3. — The extracted one-dimensional spectra showing the 
regions of detected emission lines We separate these detections to 
different knots identified in Fig. 1, mainly knot C, brightest of 
the features, and the combination of knots A+B. The spectra from 
knot C are plotted on the left while the spectra from knots A+B are 
plotted on the right. Emission lines from knot D were undetected or 
confused with the continuum emission from the foreground lensing 
galaxy in the dispersion direction of the two grism observations. 

spectra using the known redshift of 2 = 1.027 for H1429- 
0028. We made use of custom IDL scripts as well as 
the Pyraf task SPLOT with ICFIT to extract the emis¬ 
sion line flux densities and their errors from the one¬ 
dimensional spectra. To account for the contamination 
from the foreground lensing galaxy, we also model a con¬ 
tinuum and account for the contamination in our line flux 
densities. Detailed lens modeling of H1429-0028 in four 
bands have shown clear evidence for differential magnifi¬ 
cation. For the rest-frame optical observations as is the 
case for onr data the appropriate magnification is 7.9±0.8 
(jMessias et all 120141 ). 

3. RESULTS 


TABLE 1 
Emission Lines 


Line 

Component 

Flux 1 

Eq. Width 2 (A) 

Ha+[NII] 

(A+B) 

45.8 ± 2.3 

108.4 ± 5.4 

Ht+NII] 

(C) 

32.9 ± 2.6 

52.4 ± 4.2 

H/3 

(A+B) 

4.5 ± 2.6 

15.1 ± 8.8 

H/3 

(C) 

1.4 ± 3.8 

2.9 ± 8.0 

[OIII]A(5007) 

(A+B) 

6.4 ± 1.8 

18.8 ± 5.2 

[OIII]A(4959) 

(A+B) 

4.1 ± 1.9 

12.6 ± 5.9 

[OIII] A(5007+4959) 

(C) 

10.8 ± 2.2 

22.6 ± 4.7 

[SII] (doublet) 

(A+B) 

10.1 ± 1.8 

23.7 ± 4.2 

[SII] (doublet) 

(C) 

8.7 ± 2.1 

14.1 ± 3.4 


1 Line fluxes are in 10“ 17 erg s -1 cm -2 , not corrected for lens 
magnification. 

2 Equivalent widths should be considered as an upper limit due 
to potential systematic uncertainties in the background continuum 
model. 


In order to derive extinction corrections and line ra¬ 
tios, we measure the line intensities for detected bright 
regions of H1429-0028. In Fig. 1 we label the bright c orri- 
ponents following the scheme of lMessias et al.l ( 2014 ). In 


the G141 grism the brightest component C and (A+B) 
as well as parts of the Einstein ring had clear emission 
detections. For G102 only components (A+B) had a de¬ 
tectable [OIII] as well as H/?. In G102 component C had 
clear detection of the [OIII] lines but only an upper limit 
on the H/? line. Given the ratio of Ha to H/3 in knots 
(A+B), and the detected value of Ha in knot C, the ex¬ 
pected value of H/? in knot C is ^ 3.2 ± 1.9 x 10 -17 erg 
s~ 1 cm~ 2 . We measured the H/? line flux density to be 
1.4 dz 3.8 x 10 -17 erg s~ 1 cm~ 2 . This value falls within 
the estimated range but is not robust enough to be used 
for scientific analysis. 

From Fig. 1 the lens models suggest two components 
for H1429-0028, one that is compact and bright and a 
second that is extended. The two components have ef¬ 
fective radii of 0.18” ±0.01 and 0.03“±0.01, for the source 
responsible for the ring and for knots, respectively. It is 
clear from the lensing models that the quadruply imaged 
knots A to D are from the same source. While it has been 
suggested that H1429-0028 is a merger, due to the pres¬ 
ence of two components in the lens model, it is not clear 
from such a model if H1429-0028 is two separate galaxies 
or if the smaller component is a high star-forming region 
within a galaxy. Our spectral line data do not have the 
adequate velocity resolution but in the future this ques¬ 
tion can be addressed with an integral field unit (IFU) 
observations. In this work, for the line ratios, we only 
study the ratios of bright knots. Thus our line ratios 
capture the physical properties within the starbursting 
compact region or a compact galaxy. For the total Ha 
flux we add the flux from each component with their cor¬ 
responding magnification from the lens model. The val¬ 
ues are a magnification factor of /i ^ 27 for the compact 
component, which contributes to the bright knots, and /i 
^ 10 for the larger component which contributes to the 
ring. We scale the observed line flux to a total estimate of 
the line intensity across the galaxy based on K-band pho¬ 
tometry of the bright components and the diffuse rings. 
This correction results in a factor of 2.6±0.1 from the line 
fluxes measured for the sum of the components A±B±C 
to the galaxy as a whole assuming that the continuum 
fluxes detected for the other components and ring scale 
as the rest-frame optical magnitudes. 

To correct for the [Nil] contamination of Ha we 
make use of two independent methods to derive the ex¬ 
pected [Nil]/Ha line ratio and average them as the fi¬ 
nal va!ue_to_iise_Jiere. Tins follows the approach given 
in Dominguez et al . (|2013[ ). The first method from 
iSobral et al.1 (j2012f T estimates the [Nil ]/Ha ratio using 
the Ha ± [Nil] equivalent width (EWs). Using the mea¬ 
sured EWs we estimate the ratio to be 0.27 ±0.07. How¬ 
ever it could be that we are overestimating the EWs in 
our line fitting procedure due to systematic uncertainties 
associated with the model for the continuum, especially 
since the continuum is dominated by the residual fluxes 
from the foreground lensing galaxy. Thus we also employ 
a second method, but we find consistent estimates on the 
[Nil] /Ha ratio. 

The second method from lErb et ahl ( 2006 ) relies on 
a relationship derived between stellar mass of a galaxy 
and the Ha/NII ratio of that galaxy. Here, instead of 
an independent estimate of the stellar mass, we make 
use of the SED modeling in Ma et al. (in prep) with a 
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Fig. 4.— Ha lumi nos ity vs . Far-infrared for H1429-0028 compared to the extinction corrected sample of IR luminous and sub-mm 
galaxies in ITakata et alJ ( 20061 b For reference we also show the corresponding SFRs based on IR luminosity and Ha luminosity to the 
top and left of the plot, respectively. We show H1429-0028 for two cases with and without extinction correction of Ha luminosity. The 
dot-dashed line represents the case that SFRs from Ha and far-infrared are equal. H1429-0028 falls below this trend line but the difference 
between IR and Ha-based SFRs is fully consistent with the observed scatter of previous measurements. 



Fig. 5.— Left: Balmer decrement vs. IR luminosity. The bac kground data are from Takata et al. (2006). Middle: Balmer decrement 
vs. Ha luminosity. The stars come from [Dominguez et al.l (12013 ) and the diamonds are from ITakata et al.l (|2006h . The contours show the 
galaxy population traced by SDSS. Right: Balmer decrement vs. galaxy stellar mass. The contours show the galaxy population traced 
by SDSS. The diamonds correspond t o star -forming galaxies of 0.75 < z < 1.5 presented in IDomfnguez et al.l 12013 ). while the squares 
correspond to z ~ 2 from ISobral et all (12012 ). We show the expected optical attenuation Ay to the right of the right panel for corresponding 
values of H^/H^The dashed line represents the intrinsic value of the Balmer decrement. 
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Fig. 6. — Top: BPT diagram with the black line separating 
the AGN and star-forming regions iKewlev et alJ ( 20011 ). The red 
point corresponds to H1429-0028, while the black points come from 
star for m ing galaxies in the redshift range 0.75 < z < 1.5 from 
IDomm guez et al. (2013) represe nted as diamonds an d two z ~ 2 
lensed star-forming galaxies from IHainline et al. (2009) represented 
as squares. The background contour s show the galaxy population 
traced by SDSS IKewlev et alJ 120011 5 . 



Log M* [M©] 


Fig. 7. — Metallicity vs. Stellar Mass. The green band represents 
local SDSS galaxies while the blue region represe nts a second or¬ 
der fi t to SDSS extrapolated towards higher SFR (IMannucci et al.l 
1201015. The orange dashed line representing z ~ 2 galaxies from 
Isteidel et al.l (120141 b The calculated value of metallicity for H1429- 
0028 representing the 03N2 calculation using the 0[III]/H/3 to 
N[II]/Ha ratio (IPettini fe Pagell 120041) . 

stellar mass of 1.9 ± 0.02 x 10 11 M 0 derived from MAG¬ 
PHYS. The resulti ng; [ Nil]/Ha ratio is 0 .27 A 0.03. Both 
iSobral et al.l (|2012l ) and lErb et al.l (I2QQ6D make use of the 
calibration method used in IPettini fe Pagel (j2004f ). The 
two estimates are consistent with each other. Note that 
the stellar mass we use is sli ghtly lower than the value 
of ~ 3 x 10 11 M 0 quoted by IMessias et al.l (j2014f ). We 
prefer to use a revised value for the stellar mass using 
MAGPHYS models as we include new optical measure¬ 
ments of the SED and, as discussed later, we find con¬ 
sistent estimates on the extinction with SED modeling 


Fig. 8.— Star formation rate (based on the IR luminos- 
ity) vs. stellar mass. For compari s on z ~ 1 SMGs from 
Mich alowski et all (120100 ITacconi et al.l 1201015 . and IBanerii et all 
(l2onh are shown. The dashed line shows the z = 1 main sequence 

when compared to the estimates based on the Balmer 
line ratios. 

For the values we d iscuss below we take the average ra¬ 
tio of [Nil]/Ha from ISobral et al.l (|2Q12l ) and lErb et a hi 
(j2006[ ) to be 0.27±0.05. Once corrected for [Nil] we find 
Ha/H [3 = 7.49 =b 4.4. Using this ratio we calculate the 
nebula r extinction E(B — V ) following iMomcheva et ah 
( 2013 ) and find it to be 0.82 db 0.50. Using iCalzetti 
(j2001f ). the corresponding optical depth ry is 2.93 db 1.9. 
The extinction is lower than the ry value of ^ 11. 2^ \ 
for H1429-0028 in IMessias et al.l (120141 ). based on the 
broad-based SED model fitting using Magphys. A re¬ 
vised model fit to H1429-0028 using new estimates of 
the background galaxy, including deblended IRAC data, 
finds Ty ^ 4.2 zb 0.4 consistent with the estimate of ry 
from the Balmer line ratios (Ma et al. in prep). 

In Fig. 4 we compare the extinction-corrected Ha lu¬ 
minosity of H1429-0028 with other star-forming galax¬ 
ies. All data points from the literature (following Takata 
et al. 2006) are extinction corrected though we show 
the case for H1429-0028 with and without extinction 
correction. Though the apparent luminosity of H1429- 
0028 corresponds to that of a hyper-luminous infrared 
galaxy with Lir rsj 10 13 Lq, the intrinsic luminosity, 
once corrected for lensing magnification, is that of a 
ultra-luminous infrared galaxy (ULIRG). The galaxy 
falls between th e SMG s and local ULIRG s stud ied by 
iSwinbank et al.l (|2004f ) and iTakata et al.l (|2006l ) with 
rest-frame optical spectroscopy at Keck and Subaru, 
respectively. We find the extinction-corrected SFR of 
H1429-0028, at 60 d= 50 M© yr -1 , to be lower than the 
instantaneous SFR implied by the total IR luminosity, 
with a value of 390 =b 90 M© yr -1 using the iKennicuttl 
(jl998[ ) relation. Given the scatter observed in Fig. 4, 
however, we do not find this difference to be statistically 
significant. 

In Fig. 5 we compare the Balmer decrement of H1429- 
0028 against Ha, IR luminosity, and stellar mass for a 
sample of galaxies. As shown in Fig. 5 middle panel, 
for the sample of galaxies with both H /3 and Ha mea- 
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surements in the literature, we find that the extinction- 
corrected Ha lu minosity of H142 9 -0028 to be among 
the highest. The lDommguez et al.l (|2013l) sample comes 
from HST/WFC3 grism observations of z ~ 0.75-1.5 
galaxies. The SDSS-detected star-forming galaxies have 
Lhc* < 10 42 ergs s _1 , while for H1429-0028 Lh« > 10 43 


.-i 


ergs s _ \ This is consistent with the fact that H1429- 
0028 is an ULIRG. The right panel shows the trend in the 
Baimer decrement with the stellar mass such that there 
is a slight decrease in the H/? to Ha ratio with an increase 
in the stellar mass. T he plotted points a re the sample- 
averaged values from iDommguez et al.l (j2013f ) as dia¬ 
monds and lSobral et alT ( 20121) as squares in three stellar 
mass bins in both studies. These data mainly probe the 
stellar mass below a few times 10 10 M 0 . H1429-0028 is 
massive with M* ^ 10 11 M 0 and has a Balmer decre¬ 
ment that is lower than the typical star-forming galaxies 
in the same redshift range of 0.75 to 1.5. 

In Fig. 6 we compare the line ratios of [OIIIJ/H/? vs. 
[SII]/Ha. This is a variant of the more traditional BPT 
diagram (B aldwin et al.l fl~981h that involves [OIIIJ/H/3 
vs. [Nil]/Ha. Given that [Nil] is blended with Ha in 
our low-resolution data we make use of [SII]/Ha ratio. 
We make extinction corrections for the [SII]/Ha ratio 
here given the two [SII] lines are somewhat separated in 
wavelength from Ha. In this diagram H1429-0028 is con¬ 
sistent with the low metallicity end of the star forming 
regions although the ratios have large uncertainties asso¬ 
ciated with measurement errors. The measurements are 
incompatible with AGN regions of galaxies from SDSS 
data at z < 0.3. While H1429-0028 is luminous this is 
primarily due to gravitational lensing; the intrinsic lu¬ 
minosity of H1429-0028 is compatible with a galaxy star 
forming at a rate of 200 to 400 M 0 yr _1 . The lens mod- 


are compatible with 
a merger system. Interestingly a value for [Nil] /Ha of 
0.27 =b 0.03 is higher than the average [Nil ]/Ha ratio of 
0.19 d= 0.05 for the ga laxies classified as star -forming in 
the SCUBA sample of lSwinbank et all (|2004fl . and lower 
than the average for SMGs hosting AGN of 0.41 d= 0.05 
from the same study. 

In Fig. 7 we make use of the nebular line ratios, 
with the estimate of [Nil ]/Ha ratio, to make an esti¬ 
mate of the metallicity. Instead of an estimate based 
on [Nil ] /Ha ratio alon e, we make use of the 03N2 ratio 
([Pettini &; Pagell l2004f ) as the estimator here as it also 
involves the measured [OIII]/H/3 ratio. The metallicity 
value, as measured in terms of 12 ± log (O/H) was found 
to be 8.49±0.16. In Fig. 7 we compare the metallicity vs. 
the stellar mass. The figure shows the av erage metallicity 
vs. st ellar mass relations fo r both local (Mannucci et al.l 
l2010h and z ~ 2 galaxies ([Steidel et alll2014} ). H1429- 
0028 has a metallicity comparable to galaxies at z ~ 2 
despite being at 2 = 1.027. H1429-0028 has a high 

SFR, but shows no indication that it is hosting an AGN. 
H1429-0028 is metal poor despite its high stellar mass 
and argues for a scenario that it is still under a rapid 
formation phase. 

Finally in Fig. 8 we show the location of H1429-0028 
in comparison to the main sequence of galaxies at z ~ 


1 . Here we plot the total IR luminosity-based SFR of 
H1429-0028 vs. stellar mass. We fi nd that H1429-0028 
is above the z ml correlation from lElbaz et al ] (l2007h . 
For comparison, we also show other dusty star-forming 
galaxies at z ~ 1 from the literature. Fina lly the average 
gas fr action Mism /(M 0 ±Mjsm) for the iTacconi et al.l 
( 2010 ) sample of star-forming galaxies is 0.34%. The 
gas fraction for H1429-0028 is 0.25 ± 0.1%, where we 
make use of the gas mass of M ism =4.6 ±1.7 x 10 lo M w 
from ALMA CO observations reported in iMessias et al.1 

(l20Tl . 

H1429-0028 is one example of a grism observation with 
HST based on a galaxy that was first selected with the 
Herschel catalog as a lensed back ground source. Based 
on lensing models (jWardlow et al.1 l2013h , we find that 
there should be roughly 0.25 deg -2 lensed starburst 
galaxies in the redshift interval of 1 to 2. In the future 
such galaxies will be automatically included as part of 
the surveys that will be done with slitless grisms on Eu¬ 
clid and WFIRST. In the 2000 deg 2 High Latitude Deep 
survey we expect WFIRST will detect close to 500 lensed 
starbursts at z ~ 1 to 3. The study we have presented 
for one lensed galaxy can then be expanded to a large 
enough sample for detailed statistical study that probes 
the internal structure of lensed starbursts. 

4. SUMMARY 

We observed the Herschel -selected gravitationally- 
lensed st arburst galax y HATLASJ1429-0028, studied in 
detail in IMessias et al l (120141) with some initial descrip¬ 
tion in lCalanog et al 1 720141 ). We present Hubble/WFC3 
G101 and G412 grisms of HATLASJ1429-0028. The ob¬ 
servations covered the wavelength regime of 0.8 to 1.7 
fi m and resulted in detections of Ha±[NII], H/?, [SII], 
and [OIII] for several bright regions of the background 
galaxy. The Balmer line ratio Ha/H/3 of 7.5 ±4.4, when 
corrected for [Nil], results in an extinction for the star- 
burst galaxy of E(B — V) = 0.8 ± 0.5. The Ha based 
star formation rate, when corrected for extinction, is at 
the level of 60 ± 50 M 0 yr -1 , lower than the star for¬ 
mation rate of 390 ± 90 M 0 yr -1 from the total IR 
luminosity. HATLASJ1429-0028 also has a low metal¬ 
licity despite its high stellar mass at the level of 10 11 
M 0 . The combination of high stellar mass, lack of AGN 
indicators, low metallicity, and the high star-formation 
rate of HATLASJ1429-0028 suggests that this galaxy is 
still going through a rapid formation. 
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